Phase separation frustrated by the long range Coulomb interaction II: Applications 
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The theory of first order density-driven phase transitions with frustration due to the long range 
Coulomb (LRC) interaction develop on paper I of this series is applied to the following physical 
systems: i) the low density electron gas ii) electronic phase separation in the low density three 
dimensional t — J model iii) in the manganites near the charge ordered phase. We work in the ap- 
proximation that the density within each phase is uniform and we assume that the system separates 
in spherical drops of one phase hosted by the other phase with the distance between drops and the 
drop radius much larger than the interparticle distance. For i) we study a well known apparent 
instability related to a negative compressibility at low densities. We show that this does not lead 
to macroscopic drop formation as one could expect naively and the system is stable from this point 
of view. For ii) we find that the LRC interaction significantly modifies the phase diagram favoring 
uniform phases and mixed states of antiferromagnetic (AF) regions surrounded by metallic regions 
over AF regions surrounded by empty space. For iii) we show that the dependence of local densities 
of the phases on the overall density found in paper I gives a non-monotonous behavior of the Curie 
temperature on doping in agreement with experiments. 

Pacs Numbers: 64.75.1g 71.10.Hf 71.10.Ca 75.30.Vn 
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I. INTRODUCTION 

In the last decades continuous progress in the charac- 
terization and preparation of complex compounds have 
produced a variety of systems with very reach phase dia- 
grams when the concentration of some dopant is varied. 
Notable examples includes doped cuprates and mangan- 
ites where one finds different phases as the electronic 
density is varied. Quite generally and in analogy with 
familiar first order phase transitions, like the ice-liquid 
transition, it is natural to ask under what conditions one 
can find ranges of global electronic density with phase 
separation among the "ta.axrf electronic phases that one 
finds in these materials .Era 

This pmblem has arisen naturally forj-jdoped Mott 
insulatoraj'cl and Fermi liquid instabilitiestl in the con- 
text of the cuprates and also in_the related problem of 
doped magnetic semiconductors. eI 

It is by now settled that close to the Mott transition 
there is a natural tendency for the system to phase sepa- 
rate in insulator and metallic phases with different densi- 
ties. This tendency is frustrated by the long range inter- 
action which tends to favor uniform phases. As a result 
the system can choose to phase separate in a inhomoge- 
neous state with islands of one phase in the other phase 
keeping long scale neutrality. The same phenomena can 
occur in a variety of situations and in particular in the 
doped magnetoresistant manganites phase separation at 
various |Scales is observed in different regions of the phase 
digramfl Also there is evidence that the 2-dimen^onal 
(2D) electron gas phase separates at low densities.Q 

In paper I of this series (hereafter referred as I) we 



have presented a theory of phase separation frustrated 
by the long range Coulomb (LRC) and in the presence 
of a surface energy cost. We showed that if the Coulomb 
and surface energy cost are not two strong the phase sep- 
arated state survives but it is inhomogeneous. In certain 
global density range drops of one phase [A) are formed 
and hosted by the other phase {B). The free energy per 
unit volume reads: 



/ (1 - x)fA{nA) + xfEins) + er. 



(1) 



where x is the volume fraction. The first two terms are 
the bulk contribution of the A and B phases and the last 
term is the mixing energy 



1/3 



(2) 



where u{x) is a geometric factor which in the case of 
drops is: 



,(x) = 3^/3 (^) x(2 - 3x1/3 + xf/^ 



(3) 



The mixing energy includes the surface energy cost and 
the electrostatic cost. 

In our computations we have assumed that the scales of 
the inhomogeneities is much larger than the interparticle 
distance. This study is complementary to others which 
have considered the opposite limit (frustrated phase sep- 
aration at a scale comparable to the interparticle dis- 
tance) to, explain phenomena like the striped states in 
cuprates .B'EI 
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We have considered spherical drops as done by Nagaev 
and collaborators in the context of doped magnetic seijai-i 
conductors in general and of manganites in particular B-El 
However we obtain similar results in other geometries 
like a periodic array of layers and we believe that for any 
reasonable geometry similar behavior for thermodynamic 
quantities would be obtained. 

To illustrate the generality of theses ideas in this paper 
we consider some relevant applications to open problems 
in condensed matter. 

It is well known that the-Jow density electron gas has 
a negative compressibility.ll3 We discuss the fundamen- 
tal problem of the stability of the electron gas and of 
the Wiener crystal at low density against a bubble phase 
(Sec. H). The system is shown to be stable against this 
kind of phenomena showing explicitly that a negative 
compressibility can be observed in this system because 
the LRC interactions makes it stable. Interestingly nega- 
tive-eompressibility has been measured in the 2D electron 
gas. til 

To make a link with the problem of phase separation in 
doped Mott insulators wc consider frustrated PS in the 
t — J model (Sec. HI). This is one of the simplest models 
used in the context of high temperature superconductors 
where frustrated phase separation is believe to play an 
important role. We illustrate the importance of the LRC 
forces in determining the phase diagram. 

Finally we study the problem of the phase separation in 
the Manganites between a ferromagnetic metallic phase 
and a charge ordered phase (Sec. IV). This problem il- 
lustrates nicely how the anomalous behavior of local den- 
sities found in I can reflect in measurable quantities. We 
propose an explanation for the anomalous behavior of 
the Curie temperature close to a charge ordered instabil- 
ity, i.e. the Curie temperature decreases as the charge 
ordered instability is approached. We conclude with a 
summary of the main results (Sec. 0). 



II. STABILITY OF THE JELLIUM MODEL 

Here we discuss the case of a system of electrons in 
a uniform rigid background usually called the "jellium" 
model. Although we find that drops do not form in this 
case, this first discussion is very useful to illustrate the 
range of applicability of the present ideas. 

The problem is the following: It is well known that 
a low density |-electron gas has a negative electronic 
compressibility.ll3 Will this lead to drop formation?. 

To describe in an approximate way the electronic en- 
ergy one can use the Wigner interpolation formula for 
the correlation energy. In this approximation the ground 
sta±n energy per particle at zero temperature is given 

byJia 



Eel/Ry = 



2.2099 0.9163 



0.88 



where the first term is the kinetic energy, the second 
term is the exchange energy, and the last term is the 
correlation energy. Here Vg = [3/(47rn)]^/'^/ao and oq is 
Bohr's radius. 

The energy per unit volume is fei{n) — EeiJi. In Fig. ^ 
we plot fei and Eei as a function of density. These curves 
can be interpreted in two different ways. If the back- 
ground compressibility is given only by the electrostatic 
self energy (already included) then the curves represent 
the total energy of the system, background plus elec- 
trons. We call this the compressible background case. 
Two different criteria give thermodynamic instability for 
the compressible background case. First for rtog < 0.0015 
(up arrow) the compressibility is negative. More impor- 
tantly for uuq < 0.003 = Tigj (down arrow) the pressure 
is negative. The latter means that if the system is pre- 
pared with a density lower than n^j, then electrons and 
background will relax to a self bound system with a lower 
volume and n = (from now on we shall measure den- 
sities in units of Bohr inverse volume Og). We can con- 
sider this result as due to the usual Maxwell construction 
(MC) argument applied to phase separation between the 
electronic system-|-background and the vacuum. In fact 



it is easy to see that n 



satisfy a MC in which the 



MC line intersects the origin. By putting an external 
pressure densities higher than nj^j become physically ac- 
cessible. The numerical value of n^^ can change when 
more accurate forms of the correlation energy are con- 
sider but the qualitative picture will remain the same. 




as 



(4) 



FIG. 1. Energy per particle (left axis, solid line) and per 
unit volume (right axis, dashed line) as a function of density 
obtained from Wigner interpolation formula. The vertical ar- 
rows indicate the density at which the pressure is zero (lower 
one) and the density at which the jellium model contribution 
to the compressibility becomes zero (upper one). In the later 
case the corresponding change of curvature is almost indistin- 
guishable to the bare eye. 

The above interpretation is not useful in real situa- 
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tions where forces other than the electrostatic one can 
constrain the background to have a certain density. This 
leads to a second interpretation of the curves in Fig. ||. 
Since the background has an additional non-electrostatic 
contributions to the compressibility (for example coming 
from core-core repulsion of the atoms), the plots of Fig. |l| 
are not any more the total energy of the (neutral) system 
as function of density. 

A simple hypothesis to describe such a system is to 
assume in the model that the substrate is completely 
rigid. We call this the incompressible background case. 
The curves in Fig. ^ represent then the energy of many 
different realization of this system each with a different 
density. The energy including the electrostatic cost to 
change the density of the background from one system 
to the other but excludes the (infinite) non-electrostatic 
energy cost to change the density of the background. 

Now the total density is fixed at some value, the com- 
pressibility of the whole system is infinite and the above 
instability criteria do not apply any more. This however 
does not guarantee the stability of the system. One can 
imagine that the system may be unstable towards an in- 
homogeneous phase with electron rich and electron poor 
bubbles in the uniform fixed background. 

We analyze below the case in which the electron poor 
regions have zero electron density. In principle we can 
work as in I with a quadratic expansion of the free energy 
around the MC case however the free energy has now a 
simple form which can be dealt with analytically. 



A. Electron drops in background 

The compressible background case suggests that the 
system has the tendency to separate in electron rich re- 
gions and regions of zero electron density. We will take 
the A and B phases of I to be the background with no 
electrons and the background with an undercompensated 
density of electrons respectively. Consider first the case 
of low densities for the B phase. We can take for the 
bulk drop free energy the energy of a classical Wigner 
crystal i.e. the leading l/r^ term in Eq. (|^) for large Tj. 
This microscopic Wigner crystal should not be confused 
with the mesoscopic Wigner crystal that the drops would 
form. 

It is instructive to write the free energy in the following 
way: 



\ 1/3 

2-3(— + — 

nei J riei 



at low density) and the second term is the mixing en- 
ergy contribution computed in I. The latter contains the 
electrostatic bubble contribution and the surface energy 
contribution. The radius is not a free parameter but 
Rd = Rd{n,nei, (t) is the drop radius that minimizes the 
free energy. 

Notice that the drop of electrons is not neutral since 
the density of electrons is larger than the compensating 
background i.e. within the drops we are dealing with a 
charged Wigner crystal of electrons in contrast with the 
usual neutral Wigner crystal of electrons. On the other 
hand in the computation of the charging energies of the 
drops in I we have assumed for simplicity that the density 
is uniform inside the drop. In the appendix A of paper I 
we compute the correction to the electrostatic energy due 
to the non uniform electronic density at the microscopic 
scale as it should be for a charged Wigner crystal and 
conclude that this only changes numerical factors, which 
are not important for the present analysis. 

Since the drop radius has already been minimized one 
is left with the density inside the drops (or equivalently 
with Tej) to be minimized. Rd depends on density ex- 
plicitly and indirectly on a(riei). If the last term in the 
curly brackets grows with Uei faster than Ueir'^i-, the min- 
imum occurs at Uei n i.e. x = I. This corresponds 
to the uniform case. If instead rieir'^i grows faster one 
finds a solution with Uei — > oo, the term in the square 
brackets become a constant and clearly r^i > Rd- In 
this case the mesoscopic bubble model is clearly not ad- 
equate. In order for both term to balance exactly one 
finds that the surface energy has to fulfill the relation 
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(5) 



/vei . If one estimates the surface energy as a 
characteristic energy density {ueie^ /rei) times a charac- 
teristic length (re; ~ ^e/^^"^) '-'^'^ '^^^ conclude on dimen- 
sional arguments that this surface energy is precisely the 
one of a Wigner crystal. Smaller surface energies give 
drops which are too small for the mesoscopic treatment 
and larger surface energies give no drops at all. 

What about the other contributions to the bulk free 
energy in Eq. (^) which will become important as the 
density inside the drop becomes large? They only make 
the drop bulk term less negative so an even smaller value 
of the drop radius is needed to stabilize the drop solution. 
From this point of view we can conclude that mesoscopic 
or macroscopic drops of electron gas are not possible. 

The only dubious case could be close to a; = 1 {nei — > 
n). Since in this case the term in the square brackets 
can be very small. In principle this allows for large drops 
without paying too much mixing energy. However in this 
region Eq. (^ is not strictly valid since the volume frac- 
tion is close to 1 and one has to revert the geometry as 
done below. 



The volume fraction of the B phase x, has been elim- 
inated by using the constraint in the density given by 
n — xUei where riei = 3/(47rrgj) = Ub and Ua = 0. 
The first term in the brackets comes from the classical 
Wigner crystal contribution (the leading term in /s = /e/ 
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B. Drops of empty background (voids) 

We consider the possibility of formation of drops of 
zero electronic density (voids) hosted by electron rich re- 
gions with density rie/. 

We look again to the limit of the classical Wigner crys- 
tal. Now X will represent the fraction of empty elec- 
tronic volume. The constraint in the density is given by 
n = (1 — x)nei and the free energy reads: 



(6) 



-I- (uei - n)neiRa 



X 1/3 

n \ n 
3-3 1 



riei 



riei 



We see that if n^i — > n we can get i?^ >> r^i with a 
small surface and electrostatic energy (the last term in 
the curly brackets). Using the density constraint to elim- 
inate Ue/ in favor of x we find that for small x the free 
energy behaves as rt^(— 5rQ + 6i?^)x with ro given by 
n = 3/(47rrQ). Clearly to have a minimum for small 
X > we need Rd < t'o and the model does not apply. 
The full expression for the free energy taking into account 
the electron kinetic energy gives an even smaller sloop for 
the dependence of / on x so that an even smaller drop 
radius is obtained. We could still have drops with a finite 
electronic densities in both the drop and the host phases. 
In this way one can reduce the mixing energy because it 
depends on the density difference — [Eq. (||)] . One 
could expect to find a solution close to the critical den- 
sity for zero pressure n^j. We have searched for such a 
solution assuming Eei ^ {n — n^^iY' ■ It has higher energy 
than the uniform solution. 

We can conclude that a 3D electron-jellium model is 
not unstable towards mesoscopic or macroscopic drop for- 
mation and density regions where the compressibility of 
the electron gas is negative are physically accessibly. In- 
terestingly a negative cotapressibility is actually observed 
for the 2D electron gas.til Our result stems from the fact 
that both the energy gain coming from the MC and the 
energy cost have the same electrostatic origin. Of course 
we cannot discard instabilities that can occur at a micro- 
scopic scale. 



III. FRUSTRATED PHASE SEPARATION IN 
THE 3D t - J MODEL. 

In the last few years it has become clear that many of 
the strongly correlated models used to describe high tem- 
perature superconductors exhibits PS in some regions of 
parameter spaccElij Due to the strong anisotropy of these 
materials usually 2D models are consider. In this section 
we apply the idea of a Wigner crystal of drops to PS in 
models of strongly correlated electrons on a lattice. We 
will consider, for simplicity and homogeneity with I and 



the other sections of the paper, isotropic 3D lattice mod- 
els. We expect however that the results will remain qual- 
itative valid even for 2D models. Needless to say that the 
3D models are interesting on its own right given the large 
class of strongly correlated materials where anisotropy is 
not important like doped Ceo , magnetoresistant mangan- 
ites, etc. 

Usually in strongly correlated lattice models the 
Coulomb interaction is truncated to a distance of a few 
lattice sites and often only the on-site Hubbard U term 
is kept. The underlying assumption is that in a uniform 
ground state most of the interesting physics is governed 
by the short range interactions and that the effect of the 
long range interactions can be absorbed in the Madelung 
potential trough a proper Hartree renormalization of the 
on-site energies. However in a non uniform ground state 
the long range part of the interaction has an important 
role even at the Hartree level. A simple way to take this 
into account is to maintain the relevant short range in- 
teractions (e.g. the Hubbard J7), to evaluate the energy 
of the intrinsic ^ or i? phases and to add the electro- 
static and surface contribution of the drops to the total 
free energy. This means that we are still neglecting the 
Coulomb interaction at distances larger than the lattice 
constant a, as in the usual Hubbard model but we keep 
the Coulomb repulsion for mesoscopic distances of the 
order of the inhomogeneity scale. In other words in the 
Fourier transform of the Coulomb potential, Aire^ /q^ , we 
maintain terms with wave vector q close to g = that do 
not cancel with the background and hence give a large 
contribution to the energy. 

As an example of the relevance of this approach for 
strongly correlated systems we focus on the t — J model, 
one of the more often used models in the cuprates. The 
Hamiltonian is given by: 



H 



E 

<ij>,s 



J J2 



1 



<ij> 



where c\ ^ creates an electron of spin s on the site i. rii 
and Si being the electron number and spin operators re- 
spectively. Double occupation is not allowed and sum- 
mations are extended to nearest neighbour of a 3D cubic 
lattice.|-|Xlie, large J/t limit has been studied in detail 
in twop'EHO and more dimensions based on a large d 
expansion.cil We study the limiting case of small num- 
ber of electrons (hole doping close to one). This is not 
particularly relevant for the cuprates but illustrates the 
issue of frustrated PS in a strongly correlated system. 



A. Maxwell Construction analysis 

We start by reviewing the usual MC argumentsB0 in 
the absence of LRC adapted to the 3D case. The antifer- 
romagnetic phase at half filling, hereafter the B phase, 
can be model by an incompressible phase with one elec- 
tron per site, the density is given hy ub = n% = 
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Our densities refer to real electrons, not to holes. The 
energy is: 



(7) 



where bJ is the magnetic energy per bond. From csti- 
mateS|-e|f the ground state energy in the 3D Hcisenberg 
modelll^ one finds b — 0.550. 

Two different situations are found for the PS. For very 
large J/t one finds PS between the AF phase and the 
electron vacuum (AF+V) . In this case we call the vacuum 
the A phase with /a — 0. Reducing J/t one finds PS 
between the AF phase and a dilute metal of electrons 
(AF+M). In this case we call the metal the A phase with 
energy: 



35/3 



(8) 



Here t is the hopping matrix element and we have used 
the effective mass approximation in the dispersicp, rela- 
tion of the low density limit of the t — J modeLcl'O 
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FIG. 2. Free energy normalized to the incompressible phase 
free energy without LRC and (from top to bottom close to the 
origin) the metal with 6t = 0.6|/£|, PS between the AF and 
vacuum (AF+V), the metal with 6t — 2.1|/a|, MC phase sep- 
aration between the AF and the metal (AF+M). {n% = a = 1) 

In Fig. H we show the free energies in the absence of 
LRC and different values of t/{—fg). The total free en- 
ergy is given by Eq. (|l|) with Cm — 0. Since densities 
are assumed to be low we can neglect the short range 
interaction between the electrons. We define the number 
of electrons per unit cell v = n/mPg. In addition we set 
a = 1 and restore it when convenient for clarity. 

For very large J/t the AF+V solution is the lowest 
in energy. Indeed in Fig. ^ the upper curve correspond- 
ing to the uniform metallic energy does not intersects 
the AF + V line, / — vf% corresponding to MC be- 
tween the V — Q and v — \ points. Decreasing the value 



of J/t when the chemical potential of the metal fulfills 
A*a(0)(= —&t) < fg i.e. t > 6 J/2 the metallic free energy 
intersects the AF+V free energy at some finite density 
and the lowest energy state is achieved by doing MC be- 
tween the antiferromagnet and the metal. In this case as 
shown in Fig. H, one finds a pure metallic phase at small 
density and MC phase separation between the AF and 
the metal for larger density. In Fig. || we show the phase 
diagram deduced form this analysis. The dilutfi-metal can 
be unstable towards a gas of bound particlesD'Lj Here we 
do not consider this effect for simplicity. 

In the next two subsections we analyze the effect of the 
LRC interaction on the AF+V PS and the AF+M PS. 
Since the electronic free energy has a simple form we solve 
the equations exactly rather than making a linearization 
as in I. 



B. Drops of an incompressible phase in vacuum 
(AF+V) 

As shown above, in the absence of LRC interaction, 
this case is realized in the large J/t limit. Now we gen- 
eralize the above discussion with the inclusion of LRC 
and surface energy effects. The A phase is the electron 
vacuum (V) {ua — O and /a = 0) and the B phase is 
the AF with one electron per site, ng = = 1/a^ and 
energy given by Eq. (|^). The total free energy is given 
by Eq. (y) with the above replacements. An expansion of 
the densities around the MC solution (Sec. Ill of paper 
I) gives a trivial result since A = O (notice that k„i — O) 
and the densities are fixed at the MC values [Eq. (29) of 
paper I]. However this is a peculiar limit. In fact as we 
will show below the total free energy does not coincide 
with the MC free energy because of the mixing energy. 
Since the densities are fixed only the radius has to be 
determined which is given by Eq. (8) of paper I. 

The surface energy of the AF is given by the en- 
ergy cost to cut a bond divided the associated surface 
a = bJ /a? and the volume fraction is determined by the 
constraint x = n/n'g = i/. Inserting this in Eq. (8) of 
paper I we obtain: 



Rd 



15 



bJe-o 



-,1/3 



47r (2-3i^i/3 + j.)e2/a 



(9) 



As long as the dielectric constant is sufficiently large the 
radius is much larger than the lattice spacing and our 
approximations are valid. 

By writing the free energy in dimensionless variables 
we can define a coupling constant that will determine 
the transition form the AF+V solution and the AF+M 
solution in the presence of LRC. It is given by: 



a 



- fO 

JB 



97reV^« )^ 
5eo 
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(10) 



find: 



Inserting the parameters for AF drops in Eq. ( px| ) 

,rl- 



we 
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i.e. the ratio of a Coulomb energy to a magnetic energy. 

Imposing that Rd > a for = one finds a < 3/2 so 
we will concentrate on this range of coupling. 

From Eqs. (|l|),(|^) we obtain the free energy as a func- 
tion of density: 

= [-^+ (24^)'^'""MK-/b) (12) 

The first term in the square brackets is the bulk contri- 
bution and the second term is the mixing energy con- 
tribution. This is plotted in Fig. ^ for different values 
of the coupling a. As in the above the results are only 
rigorously valid for small x { ^ v ). 
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V 

FIG. 3. f/{-f%) for a phase of drops of an incompressible 
phase in vacuum for a = 0, 0.5, 1, 2 (from bottom to top) vs. 

V. 

The effect of the mixing energy is to bend upwards the 
AF-I-V free energy (see Fig. ||) so that the metallic phase 
can become stable with a lower value of t with respect 
with the case with no LRC force (compare with Fig. |^). 
One can show that for a > 1 the AF-f V solution is never 
stable and one has either a uniform metal of an AF-I-M 
solution depending on doping. For a < 1 the drops can 
coexist with a metal or not depending on the value of 
t/J. We will analyze the competition with the AF-I-M 
solution in the next section. 



C. Drops of an incompressible phase in a metallic 
host (AF+M) 

Reducing J jt at some point the solution of the previous 
section (AF-t-V) will not be stable any more. This has 



already been shown in the absence of LRC interaction 
(a = 0). We consider now the A phase to be the metal. 

The surface energy will have now density dependent 
contributions coming from the metal. However, since we 
are in the low density limit the surface energy will be 
dominated by the magnetic surface energy described in 
the previous case and can be taken as constant. The 
AF-f-V solution is then not any more stable when (0) = 

-6i<MAF+v(0) = («-l)(-/l^). 

In Fig. H we show the locus of stability of the AF-I-V 
solution in the t — a plane. Above the line the stable 
solution is either a uniform metal or drops of AF in the 
metal depending on the density. In Ref. the ratio 
of J jt below which the AF-I-V solution is not stable for 
a = is called Y^. In 2D they found Yc{^) = 3.4367 
and yc(O) 4 for d — > oo.EJ Using their estimate of 
the 3D AF energy we have yc(0) = 3.637. Fig. ^ shows 
that Yc (proportional to the critical value of —f%/&t) 
increases with a. Remarkably in the presence of LRC 
forces a smaller t is enough to stabilize metallic phase 
regions. In other words we can have a situation in which 
without LRC forces all the electrons are in a self bounded 
AF state and as the LRC forces are switched-on some 
electrons "evaporate" to form a dilute gas around the 
AF drops. 



1.2 




a 



FIG. 4. Locus of stability of the AF drops in vacuum 
(AF-I-V) in the t — a plane. Above the line the more sta- 
ble solution depends on density. 

To solve for the AF-I-M drop solution the free energy 
now has to be minimized with respect to the radius and 
the density of the metal subject to the constraint n — 
x'nPg -|- (1 — x)nA- We are implicitly assuming that the 
density is not low enough to form a Wigner crystal of 
electrons. One can check that for reasonable parameters 
and increasing a the radius of the drop becomes of the 
order of the lattice constant much before an electronic 
Wigner crystal can form. 

Above the boundary line on Fig. one finds either 
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a uniform metal or AF+M depending on density. This 
can be seen in Fig. ^ where we plot the free energies 
for 6t = 2.1|/g| and a = 0.5. Above a certain density 
"nbi f = /a^ we have coexistence of AF drops in the 
metallic host. The behavior close to ribif is very similar 
to the behavior for parabolic free energies in the UDA of 
paper I. 



As the B density grows the nA density decreases due 
to the effect discussed in paper I. The B density is kept 
constant at due to the incompressibility. In Fig. Q 
we show this behavior. In real systems this effect can 
be detected through physical properties which depend on 
the local densities as is disused below for the manganites. 
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FIG. 5. Free energy normalized to the incompressible phase 
free energy with parameters a = 0.5 and 6t = 2.1|/^|. We 
show the metallic free energy, the AF+V free energy and the 
AF+M free energy. The cross indicates the value with a; = 
of the AF+M drop solution. 

Here also there is a bifurcation of the solution and in- 
creasing the density, the AF drops appear with a nonzero 
value of the volume fraction. However with the present 
parameters the initial volume fraction is very small (See 
Fig. I). 



FIG. 7. Density in the metallic host vs. total density u for 
&t — 2.1 1 /a I and different values of a. 

In Fig. H we show the phase diagram in the absence 
of LRC force {a = 0) and for a — 0.5. We see that a 
portion of the phase diagram in which a uniform solution 
is unstable towards PS without LRC, for a > becomes 
stable and the AF+M solution extends its region with 
respect to to the AF+V solution due to the "evaporation" 
effect. 




FIG. 6. Volume fraction vs. v for AF drops in a metallic 
host and parameters as in Fig. |^. 
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FIG. 8. Phase diagram for the large J/t hmit of the 3D 
t — J model without LRC (q — 0) and with a small LRC 
{a — 0.5). The high density part has to be taken with care 
since the density of electrons can be large in the metal so that 
short range interactions within the metallic phase cannot be 
neglected any more, and also drops of AF loose sense. 

IV. APPLICATION TO THE MANGANITES 

As a further application we consider the magnetoresis- 
tant manganitesEl hke Lai_j,Caj,Mn03. In the last years 
strong experimental evidence has accumulated indicat- 
ing that inhomogeneous phase separation occursjpj-these 
materials in certain regions of parameter space.li3"EEl 

At y = all Mn have formal valence 3+. Each ion has 
three electrons in t2g orbitals and one electron in an Cg 
orbital. The four electrons spins are all parallel due to 
the strong Hund's rule coupling forming an = 2 spin. 
The spins of different Mn^+ ions form an antiferromag- 
netic (AF) phase due to the superexchange interaction 
and the system is an insulator. As the eg band is doped, 
mobile holes tend to align the S = 3/2 core (t2g) spins in 
different ions because this maximizes the transfer integral 
and minimize the holes kinetic energy leading to a fer- 
romagnetic state.p-This is the so called double exchange 
(DE) mechanism.Ero 

Experimentally one finds indeed the ferromagnetic 
state but close to y = 0.5 a new charge ordered (CO) 
insulating phase with a chessboard structure of Mn'^+ 
and Mn^+ is stabilized. The CO phase is not predicted 
by the conventional DE but does appear in more reeent 
theories incorporating Mn^Ain Coulomb repulsionc3 or 
orbital degrees of freedom.o 

Close to y — and y = 0.5 the metallic ferromagnetic 
(FM) phase competes with the corresponding insulating 
phase. The drop state due to the competition between 
the y = AF phase and the metallic phase taking into 
account the LRC interaction has lW(3jn studied theoreti- 
cally by Nagaev and collaborators ffl Evidence for such 
a phase has recently been found in neutron scattering 
experiments.113 

Here we will analyze the competition of the CO phase 
with the FM phase close to y — 0.5 and show that a 
phase separated state can explain the puzzli ng . m aximum 
of the Curie temperature Tc{y) at y ^ 0.35.l3O On the 
contrary, conventional DE would predict that the Curie 
temperature is maximum at half doping {y — 0.5) be- 
cause for this filling the kinetic energy of the holes is 
maximized. 

We will consider a mixed state in which the A phase 
is the ferromagnetic metal (FM) and the B phase is the 
charge order state at y — 0.5 which corresponds to in- 
verse specific volume nP^ — 0.5/a^. In the following the 
densities refer to holes (i.e. the concentration of Mn'*+ 
ions). 

In the FM phase the core spins of the Mn ions are 
fully polarized and the mobile holes have the maximum 



bandwidth W. In ordec.to model the FM in a simple 
fashion we follow VarmaEa and take a flat density of states 
with bandwidth W. The FM free energy at T = is then 
given by the cohesive energy of the holes in the fully 
polarized state: 

fA{nA) = ^{nA-nlf. (13) 

We have chosen to measure the single particle energies 
from the center of the band and we have dropped a con- 
stant which can be absorbed in the free energy constant 
of the B phase . At finite temperatures one has to con- 
sider the entropy contribution to the free energy. How- 
ever for a given temperature one can expand the full A 
free energy around the n% density and an expression like 
Eq. (|l^ is still valid with an effective temperature de- 
pendent W . 

The CO state can be modeled as a doped incompress- 
ible phase around the inverse specific volume n^. The 
free energy at T = can be put as: 

fB{nB) ^^Wb~ n%\ + eoiuB ~ n%) + f%. (14) 

measures the difference in free energy per hole be- 
tween the CO state and the FM state a.t y = 0.5 {n = n^) 
and eo controls the difference in chemical potentials of the 
two phases. Eq is the gap in the charge order state, i.e. 
the difference between the energies to create defects with 
one added hole and one removed hole without destroying 
the CO state. (It should be of the order of the activation 
energy in the transport properties of a pure CO state). 
The dip in the free energy at n = will become rounded 
with temperature. For simplicity we will neglect this ef- 
fect. For temperatures much smaller than the gap this 
is a good approximation. Even if the temperature gets 
comparable to the gap a small rounding of the CO free 
energy close to ub = will not affect significantly the 
results close to the density at which drops first appear 
(ntif). 

The chemical potential of the CO state at T = is 
given by: 

H+^eo + Eg/2 nB>n% (15) 

^l-=eo-Eg/2 nB<n%. (16) 

By construction the discontinuity at is equal to the 
gap as it should be. is the energy to create a Mn'^+ 
defect in the CO state, i.e. it is the energy to remove 
a hole in the CO state. This single particle energy is 
measured from the same reference energy as the one used 
for the A phase in Eq. (p^. This fixes the value of cq. 

In Fig. H we have plotted /s as a function of y. 
We constructed the free energies for the uniform phases 
phenomenologically, by relying on the metallic and insu- 
lating character of each phase and on the fact that due 
to the different magnetic symmetry they cannot be joint 
with continuity but a level crossing should occur as a 
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function of y. It is interesting to note that a recent mi- 
croscopic model gives pijactically the same energy scheme 
as a function of doping.ca 

In Fig. H we report also the MC and the free energy for 
the drop solution for a = 3 and W' = Wn^/^—fg ) — 
15 (thick line). The coupling constant a is defined in 
Eq. ( JlCl) . W measures the effective bandwidth in units 
of (-^). 




y 

FIG. 9. Free energy normalized to the incompressible phase 
free energy at y = 0.5 with a bandwidth W — 15\fg \/n%. We 
show the FM free energy {/a), the CO free energy with a 
large negative value of the FM+CO free energy for a = 
(MC), and for a = 3. The cross indicates the value with a; = 
of the FM+CO drop solution. The arrow indicates the same 
for MC {yo). 

For the sake of simplicity we assume that the slope 
of /b for ub < rfi^ (i-'^- ^^~) large that Jb never 

crosses the drop solution. Under these simplifying con- 
dition one of the phases involved in phase separation is 
always the defect free CO state at y = 0.5. In this sit- 
uation Eg and cq do not enter into the problem and a 
precise description of fB{n-B) is not needed. For example 
we have neglected the kinetic energy of the defects which 
will give some curvature to /B(n_B) but will not change 
the present picture. 

Alternatively to a we could use the coupling constant 
A defined in Sec. Ill of paper I since the FM free energy 
is parabolic (fc^^ = = Wa'^) and the CO free energy 
can be consider as the fcs ^ limit of a parabola. The 
two coupling constants are related by 



A = 



22/3 



3 {wy/^' 



(17) 



Specifically A = 0.64 for a = 3 and W = 15. Notice 
however that here (as in the previous section) we can 
introduce two dimensionless parameters. One is a (or 
alternatively A) and the other is W' . The latter plays the 



same role as t/ J in the previous section. In particular it 
fixes the MC densities as follows. In the absence of LRC 
interaction {a = 0) and for W' < 4 the /a parabola is 
too flat and PS between the FM and the CO state is not 
possible. One gets PS between the CO state and v acuum 
(this is similar to the AF+V PS considered in Sec. Ill 
For W' > 4 Maxwell construction gives PS between FM 
and CO with the critical doping given by: 



1 



(18) 



In Fig. ^ the value of yo is indicated by an arrow. 
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FIG. 10. Volume fraction vs. y for CO drops in a FM host 
and parameters as in Fig. |^. The lower branch close to y ~ 0.4 
is unphysical. 

In the presence of LRC interaction the range of coexis- 
tence contracts with respect to the Maxwell construction 
case. The transition from the FM to the drop solution is 
quite abrupt at ytif — 0.38 (for a = 3) with a substantial 
jump of the volume fraction from zero to a finite value 
XMf = 0.17. (See Fig. |l|). 

In Fig. |ll| we show the local density inside the metal- 
lic region. For y < ytif the stable phase is uniform FM 
and the total density coincides with the nominal density 
n = y/a^. For y > yuf the drop solution is stable and the 
density in the metallic region decreases with increasing 
nominal density. As discussed in paper I, in deriving this 
effect it is important that the density dependence of the 
surface energy can be neglected. The strongest depen- 
dence of the surface energy is expected to arise from the 
kinetic energy of the metal. However, this dependence is 
important close io y ^ \ and y ^ Q and can be safely 
neglected close to y ^ 0.5. 

The decrease of local density with increase of global 
density can explain the non-monotonous dependence of 
Curie temperature on y. Since the ferromagnetic 
interaction between the core Mn spins is mediated by 
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the conduction electrons through the double exchange 
mechanism one expects the Curie temperature to be a 
monotonous increasing function of the local metallic den- 
sity n^(< n^) of the FM phase. We associate the re- 
gion in which the Curie temperature increases with dop- 
ing, i.e. the "normal" region (roughly 0.1 < y < 0.33 
for Lai-yCaj^MnOs), with a uniform FM phase and the 
"anomalous" regions in which the Curie temperature de- 
creases with doping with a drop state. In the latter state 
ua decreases with doping and this gives the anomalous 
behavior of as a function of doping. 




FIG. 11. Density in the metallic host vs. doping y for 
parameters as in Fig. ^. 

To be more specific we assume the following simple 
form for the dependence of the Curie temperature on the 
local FM density: 



ic(0.5)-ic(0) 



4(1 - a^UA) * a^UA 



(19) 



We are using uncapitalized t for the local Curie temper- 
ature of the FM phase to distinguish it from the true 
Curie temperature of the system which is a function of 
the overall doping Tc{y). A similar form as Eq. ( [l9| ) 
with tpifi) = was derived by Varma for a uniform FM 
phaseEa More sophisticate treatments also give a form 
roughly parabohc with tc{0) > O.c^ 

For a uniform FM phase y = uao,^, T^c{y) = tciu)- This 
fits correctly the experimental data in the normal region. 
We can use this fit to fix the parameters in Eq. (^9|) . For 
Lai-yCa^MnOa one obtains tc{0) ^ 80K and ic(0.5) ^ 
300K. Close to y = 0.5 we have the anomalous behavior 
and the measured Tc differs considerably from tc{0.5), the 
Curie temperature of an hypothetically uniform phase. 
For example, experimentally TjO.b) ~ 225K. 

From the known nA{y) (Fig. ^l|) and Eq. ( p^ ) we com- 
pute Tciy) = tc[a^nAiy)]- This curve (which is quite sim- 
ilar to the experimental one) is shown in Fig. O. Indeed 



we see that the drop solution combined with the uniform 
solution for y < yi,if gives a non- monotonous behavior 
of Tc{y). In evaluating the theoretical curve to be com- 
pared with the experimental data we fix the values of a 
and W in the following way. We associate the experi- 
mental maximum in with the bifurcation point, i.e. 
the doping at which the uniform solution switches to the 
drop solution. This gives us an experimental value of the 
bifurcation doping y^f ^ 0.38. From the experimental 
data we also obtain the depression of the Curie tempera- 
ture [Tc(0.5) - ic(0)]/[ic(0.5) - ic(0)] ~ 0.66. With these 
two dimensionless numbers we obtain the dimensionless 
parameters of our theory and find a = 3 and W — 15 
i.e. the values that we have been using in the present sec- 
tion. A rough microscopic estimate of these parameters 
is given in Appendix ^ to show that indeed the above 
values are reasonable for the manganites. 

A similar behavior as the one discussed here for y ~ 0.5 
is observed close to ?/ ~ 0. We speculate that this is due 
to the same general phenomena involving inhomogeneous 
phase separation between the insulator at y = and the 
ferromagnetic metal. 
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FIG. 12. Tc of the FM host minus tc(0) normalized to 
tc(0.5) — tc{Q) as a function of doping y. Parameters are as 
in Fig. ^. We show the Tc in the uniform solution and in the 
drop solution. 



V. CONCLUSIONS 

In this work we have applied the ideas developed in 
I to three different physical systems. First we analyzed 
well known apparent instabilities in the low density limit 
of the jellium model. Usual instability criteria like a neg- 
ative compressibility are formulated for a neutral system 
and should be taken with care in a charged system with a 
compensating background. In a charged system an insta- 
bility of the kind implied by a MC analysis is a necessary 
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but not sufficient condition for mesoscopic phase separa- 
tion. In fact we have seen that for the Wigner crystal of 
electrons the Coulomb strength and the surface energy 
balance in such a way that large drops are not possi- 
ble. This is basically due to the fact that the energy gain 
from MC and the energy cost due to the LRC interaction 
and surface energy have all the same electrostatic origin. 
This prevent the mixing energy to be small compared 
to the MC energy as required in I to have mesoscopic 
PS. We mention that in the 2D electron gas (not con- 
sider here) there is evidence both for a density region 
of negativej-compresibilities and at lower densities phase 
separation.ElEil For the latter however the effect of disor- 
der not consider by us may be crucial. 

The study of t — J model illustrates the importance 
of the mixing energy on determining the phase diagram. 
The LRC interaction tends to stabilize the non-separated 
uniform phases in the presence of a rigid background. 
Apart from this intuitive effect the LRC interaction can 
also favor one mixed state over another. In fact we 
showed that the LRC interaction can transform the clus- 
ters of AF electrons in vacuum into clusters of AF elec- 
trons in equilibrium with its vapor (a dilute electron gas), 
a phenomenon which we referred to as "evaporation" . 

In paper I of this series we showed that the local den- 
sities of the two phases tend to have an anomalous be- 
havior in the mixed state. Both of them tend to decrease 
when the global density increases. This non-linear effect 
can affect properties of the system which are sensitive 
to the local density, as we have illustrated for the Curie 
temperature in the manganites. We have thus provided 
an explanation to anomalies that occur in the phase dia- 
gram, i.e. A decreasing Curie temperature when the CO 
state is approached by varying the doping. 



APPENDIX A: MICROSCOPIC ESTIMATE OF 
PARAMETERS IN THE MANGANITES 



In Sec. IV we find that the parameters A = 0.64 and 
W = 15 give a curve Tc{y) similar to the experimental 
one. To decide if these parameters are reasonable one 
needs a microscopic computation. 



To evaluate W which appears in Eq. (18) we refer 
to a recent zero temperature microscopic computation 
which taites into account double exchange and orbital 
ordering.E3 Their Fig. 2 showing the free energy (with- 
out LRC) is quite similar to our a = curves in Fig. ||. 
From there we take ~ 0.24 which determines W ^ Ih 
[Eq. ([ij)] in agreement with the value we used to fit 
Tc- A is more difficult to obtain because it requires a 
microscopic computation of surface energies and screen- 
ing effects. We parameterize the surface energy by a di- 
mensionless quantity 7 defined by ct = "fW/a?' . Putting 
60 = (0.5 - yo)/a^ and k^^ = k'^ = Wa^ in Eq. (25) of 
paper I we get: 
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1/3 



v2/3 



(0.5-2/0)^/3 VoaW 



1/3 



(Al) 



For the bandwidth we can take an estimate based on 
Mattheiss's local density approximation (See Ref. [2^ ) 
W — 2.5eV. For a cubic array of Mn with a Mn-Mn 
distance of 4Awe get e^/a = SAeV for the bare Coulomb 
strength. Inserting the numerical values in the above 
equation we have 



A - 21 



2\ 1/3 



One obtains A ~ 0.64, the value we have used in Sec. IV, 
by taking eo ~ 100 and 7 ^ 0.05. These are reasonable 
parameters considering that eq should be understood as 
a static dielectric constant taking into account conven- 
tional dielectric screening plus Thomas-Fermi screening 
effects (Sec. IV of paper I) and should be a small 
fraction of the bandwidth. 

We mention that since a real background is never per- 
fectly rigid, a volume relaxation will also occur inside 
the drop phase. In general the positive background will 
contract in the electron rich phase and expand in the 
electron poor phase to reduce the mismatch between the 
ionic positive density and the electronic density. This is 
in agreement whit the situation in Pro.rCao.sMnOs where 
the electron poor CO phase has a larger volume than the 
electron rich FM phase. EZl Clearly this effect has to be 
included in the effective definition of eo- 

The drop radius reads: 



321/3 



A[(5a3)2(0.5 - yo)^{2 - 3a;i/3 + 2;)]i/3 



(A2) 



Using the above parameters we can estimate the radius 
at the onset {xbif — 0.17) to be of the order of 

Rd - 10a 

Correspondingly the cell radius is Rc = Rd/x^^^ ^ 18a. 
We see that these scales are much larger than the lattice 
constant and our approximations apply. 
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